Relativistic model for nuclear matter and atomic nuclei with momentum-dependent 

self-energies 
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The Lagrangian density of standard relativistic mean-field (RMF) models with density-dependent 
meson-nucleon coupling vertices is modified by introducing couplings of the meson fields to derivative 
nucleon densities. As a consequence, the nucleon self energies, that describe the effective in-medium 
interaction, become momentum dependent. In this approach it is possible to increase the effective 
(Landau) mass of the nucleons, that is related to the density of states at the Fermi energy, as 
compared to conventional relativistic models. At the same time the relativistic effective (Dirac) mass 
is kept small in order to obtain a realistic strength of the spin-orbit interaction. Additionally, the 
empirical Schrodinger-equivalent central optical potential from Dirac phenomenology is reasonably 
well described. A parametrization of the model is obtained by a fit to properties of doubly magic 
atomic nuclei. Results for symmetric nuclear matter, neutron matter and finite nuclei are discussed. 

PACS numbers: 21.60.-n, 21.30.Fc, 21.65.-ff, 21.10.-k 



I. INTRODUCTION 

The relativistic description of nuclear matter and 
atomic nuclei advanced considerably in recent years by 
improving the underlying density functional 
Originally, these nuclear systems were treated in the lan- 
guage of a relativistic quantum field theory (Quantum 
Hadrodynamics) based on a relativistic Lagrangian den- 
sity that contains nucleons and mesons as degrees of free- 
dom. In recent years, the viewpoint has changed and the 
approach is considered as an effective field theory, see e.g. 
P . The model can be seen as a special variant of density 
functional theory and a connection to the underlying 
theory of Quantum Chromodynamics was established [3 ■ 

In principle, it is possible to describe the ground state 
of nuclei exactly due to the Hohenberg-Kohn theorem 
as soon as the correct functional of the local baryon den- 
sity is known. Unfortunately, there is no explicit scheme 
for constructing such a functional. Nuclear theory is con- 
cerned to develop a sufficiently accurate approximation 
guided by experience and the physical knowledge of the 
system. Instead of using functionals based on the den- 
sity alone, one can start from a functional that contains 
additional densities, e.g. the kinetic energy density. This 
approach automatically leads to the Kohn-Sham equa- 
tions for the nucleons i.e. a Schrodinger-type equation 
in the non-relativistic case and a Dirac-type equation in 
the relativistic case, respectively. In these equations the 
effective interaction is expressed by potentials or by self- 
energies, respectively, that depend non-trivially on the 
nucleon fields. 

There is a big advantage of the Kohn-Sham approach 
since certain physical effects are more easily incorporated 
into the description. Functionals that are solely based 
on the density can be compared to a description of fi- 
nite nuclei in a Thomas-Fermi approximation where no 
shell effects appear. Considering also the kinetic density 
in the functional, corresponding to a Hartree or Hartree- 



Fock description, shell effects emerge automatically. In 
non-relativistic approaches a spin current density is in- 
cluded in the functional to generate the large spin-orbit 
splittings observed in finite nuclei. These three densi- 
ties in the non-relativistic framework are replaced by the 
baryon (or vector) density, the kinetic density and the 
scalar density in relativistic approaches. Pairing effects 
can be included by introducing the corresponding pairing 
densities. Additionally, in case of neutron-proton asym- 
metric systems, one has to deal with both isoscalar and 
isovector versions of these quantities. 

A general functional can be constructed in terms of 
the considered densities and their derivatives. Due to 
the Lorentz structure and the isospin degree of freedom 
there is a large variety of possible contributions that can 
appear in a relativistic density functional and it is not 
clear what possible densities and combinations are really 
relevant in the description. The density functional can be 
expanded in powers of the nucleon field, the meson fields, 
and their derivatives. The importance of the individual 
contributions can be estimated in a systematic approach 
that is guided by the principles of naive dimensional anal- 
ysis and of the naturalness of the appearing coupling con- 
stants P. ll0lllTI| . It is also possible to compare with more 
fundamental approaches based on QCD and chiral sym- 
metry. For large densities, however, this approach might 
be problematic due to the resulting polynomial depen- 
dence of the functional on the Fermi momentum or the 
density. Well-behaved rational approximations seem to 
be a promising alternative. Extrapolations to high densi- 
ties are necessary in the application of the model, e.g. in 
the description of heavy- ion collisions or neutron stars. 

The usual starting point in relativistic models is a cer- 
tain Lagrangian density instead of an equivalent energy 
functional. This Lagrangian can be treated in differ- 
ent approximations. For a practical application to nu- 
clear matter and finite nuclei one usually employs the 
mean-field approximation. In this approach the nucleon 
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self-energies in nucler matter become simple functions of 
various nucleon densities. In principle, one can go be- 
yond the mean-field approximation in a quantum field 
theoretical treatment by a systematic diagrammatic ex- 
pansion. In this case, the nucleon self-energies will have 
a much more complicated structure that takes the non- 
trivial effects of the nuclear medium into account. In 
Dirac Brueckner calculations of nuclear matter based on 
a given nucleon-nucleon interaction it is possible to ex- 
tract the corresponding nucleon self energies, see, e.g., 
m m El HI [13. They show a dependence on both 
the density of the nuclear medium and the energy and 
momentum of the nucleon. An alternative approach to 
these "ab-initio" strategies for the nuclear many-body 
problem is a phenomenological treatment of the problem 
where the mean-field approximation is retained with its 
virtues in applying the model self-consistently to nuclear 
matter and finite nuclei with reasonable effort. However, 
the basic Lagrangian has to be modified for a quantita- 
tive description of the various properties and it becomes 
mandatory to introduce new terms or to modify existing 
contributions in the Lagrangian. The parameters of the 
model have to be fitted to properties of nuclear matter 
and atomic nuclei. They cannot be derived in a simple 
manner from a more basic description. In this work the 
phenomenological strategy is followed since it is known 
from experience that this approach is very successful. A 
selection of the relevant terms and of a more generalized 
functional form that is guided by other approaches, e.g. 
Dirac Brueckner calculations, and by physical intuition 
may lead to reasonable results. 

In relativistic approaches to nuclear structure it is com- 
mon to introduce meson fields as explicit degrees of free- 
dom in the Lagrangian or density functional. In the mod- 
ern point of view, they have to be seen as auxiliary fields 
that, a priori, have not to be identified with the corre- 
sponding mesons in free space although they share the 
same Lorentz and isospin structure. In principle, these 
fields can be completely removed from the theory as, e.g . , 
in the relativistic point-coupling models 0, 0, ^3 . 
In the following, however, we will keep these meson fields 
in the model for convenience. 

A quantitative description of nuclear systems was 
achieved by considering an explicit medium dependence 
of the effective interaction. For that purpose non-linear 
self-interactions of the mesons were introduced and the 
corresponding coupling constants were fitted to prop- 
erties of nuclear matter and finite nuclei mostly close 
to the valleyof stability^ee e.g^^fl l |2 L jl il HI 
mmiiiElllllliMEl ISIm H At large 
densities, however, this approach is not really reliable 
due to the polynomial dependence on the fields and in- 
stabilities can occur. Alternatively, a density depen- 
dence of the meson- nucleon couplings was considered 

[11 111 IM IM El 111 111 El El 113. In Ref. 113 it 

has been pointed out that, in general, the couplings have 
to be treated as functionals of the nucleon fields leading 
to rearrangement contributions to the self-energies that 



are necessary for the thermodynamical consistency of the 
model. It is possible to choose functional forms (e.g. 
well-behaved rational functions) that are motivated by 
results of Dirac-Brueckner calculations of nuclear matter 
[lli .41.J . This class of models seems to be the more flex- 
ible approach with proper high density behavior. Origi- 
nally, the medium dependence was introduced only in the 
isoscalar part of the interaction. In recent years it was 
realized that it is also necessary to allow for a density 
dependence in the isovector channel in order to obtain 
a reasonable description of the neutron skin thickness of 
stable nuclei, the neutron matter equation of state and 
a reliable extrapolation to exotic nuclei 0, E3 • In all 
these standard relativistic approaches the self-energies in 
the Dirac equation for the nucleons depend nontrivially 
on the various densities. However, they are the same for 
all protons and neutrons, respectively, independent of the 
single-particle state. 

Despite the success of the relativistic approach there 
are still some deficiencies that have to be dealt with in 
further extensions of the underlying density functional. 
The size of the scalar (S) and vector (S^) self-energies 
in the interior of a finite nucleus is well determined by 
requiring that the corresponding central and spin-orbit 
potentials have the correct strength. (For alternative 
definitions of the non-relativistic potential see, e.g., Ref. 
(48.] . ) They are derived in a nonrelativistic reduction of 
the Dirac equation. The Schrodinger-equivalent central 
potential 

Fcen = - S + (E^ - Eg) (1) 

and the spin-orbit potential 

V -^S^Ri^ziA (2) 

h + m — IjQ — 1j Zmr 

of a nucleon with rcstmass m and energy E in a spherical 
nucleus are determined by the sum and the difference of 
the time component of the vector self-energy Eg and the 
scalar self-energy S, fixing both quantities. Correspond- 
ingly, in symmetric nuclear matter at saturation density 
Psat = 2/c|,/(37r^) with Fermi momentum pp = fikp the 
relativistic effective (Dirac) mass 

m* = m — E (3) 

that is independent of the nucleon momentum and the 
chemical potential 

= Eo + ^{m*f+pl (4) 

are well determined. 

In conventional RMF models the self-energies are in- 
dependent of the nucleon energy and momentum but 
they exhibit a strong density dependence. However, from 
Dirac-Brueckner calculations one would expect that the 
self-energies also depend on the energy and momentum 
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of the nucleon in the medium. The central potential 
allows to quantify this effect in both relativistic and non- 
relativistic models. Standard RMF models show a linear 
increase of the central potential. Only models with ex- 
plicit energy or momentum dependent self-energies show 
a different result. In order to compare the energy or mo- 
mentum dependence of ^ with empirical data it is com- 
mon practice to compare the model predictions with the 
optical potential in nuclear matter extracted from Dirac 
phenomenology (DP) for elastic proton-nucleus scatter- 
ing 23l ■ In this approach scattering observables are 
fitted to experimental data up to approx. 1 GeV kinetic 
energy by varying the strength of the vector and scalar 
self-energies (real and imaginary part) of the target. The 
potential shape is described in simple parametrizations 

47] or taken from microscopic descriptions ■ 
DP the self-energies depend on the proton energy. The 
derived Schrodinger-equivalent central potential, extrap- 
olated to nuclear matter, exhibits a non- linear energy 
dependence. At high energies a saturation of the op- 
tical potential with a value around 50 MeV is observed. 
(Imaginary contributions of the self-energies as extracted 
from Dirac phenomenology have an almost negligible in- 
fluence on the real part of the corresponding optical po- 
tential. They enter only through the terms quadratic 
in the self-energies in Eq. These contributions al- 
most cancel each other since they are of similar magni- 
tude.) At low kinetic energies, the optical potential is 
a nearly linear function of the energy. The self-energies 
in DP itself show an almost linear energy dependence 
for not too high energies. This observation and the ex- 
perience from Dirac-Brueckner calculations motivates to 
consider a phcnomcnological extension of standard rela- 
tivistic models that generates energy or momentum de- 
pendent self-energies already on the mean-field level. 

Besides the relativistic effective (Dirac) mass m* there 
are a number of other effective masses that have been 
introduced in the literature. For a general overview of 
definitions for effective masses in relativistic and non- 
relativistic calculations see |49 j. An important quantity 
is the effective Landau mass 

mo^=P^ (5) 

that is related to the density of states g = mcsp / {^ivKf 
in nuclear matter. The Landau mass is easily calculated 
from the dispersion relation 

(£;-S]o)' = (m -!])'+ (6) 

for a nucleon. In relativistic models without energy or 
momentum dependent self-energies one obtains 

meff = £; - So = ^{m* f+pl (7) 

at the Fermi momentum pp. The relativistic effective 
mass cannot be adjusted arbitrarily if a reasonable de- 
scription of the spin-orbit interaction in atomic nuclei is 



required. In usual relativistic models we find m* « 0.55m 
and TOcff ~ 0.62m > m* at the Fermi momentum but the 
latter value is rather small. Comparing non-relativistic 
Skyrme Hartree-Fock calculations for giant resonances in 
the random phase approximation with experimental data 
a value 0.78m for nics was extracted |5nl |. The low Lan- 
dau mass in conventional RMF models indicates that the 
level density at the Fermi energy is too small and single- 
particle levels in finite nuclei are too much spreaded, a 
well known result of detailed calculations. 

Assuming an energy dependence of the self-energies, 
the effective Landau mass is given by 

„., = ,i5-E„)(l-§;)+(„,-E)§, (8) 

It is clear that the Landau mass can be adjusted more 
freely if the model allows for energy dependent self ener- 
gies. However, one can not expect that the density and 
energy /momentum dependence of the Landau mass as 
observed, e.g., in DB calculations can be well represented. 
The result ^ is consistent with the non-relativistic form 

^^'^ = - ^) 

with an energy-dependent optical potential V taken from 
([T|) in the non-relativistic approximation. 

There are attempts to cure the problem of the small 
Landau mass and the related low density of levels in 
relativistic models. In Ref. l5l| a linear energy depen- 
dence of the self-energies was introduced heuristically in 
an energy window around the Fermi energy in order to 
increase the effective (Landau) mass and the level den- 
sity. However, the full self-consistency and the Lorentz 
invariance of the RMF model are lost. In another ap- 
proach the Landau mass mes was increased to the value 
of 0.76m that was necessary to obtain reasonable /3 decay 
half-lives Correspondingly, the Dirac effective mass 
had to be increased by hand to m* — 0.67m destroying 
the usually good description of the spin-orbit splittings 
in nuclei without additional modifications of the density 
functional. The authors were able to compensate this re- 
duction of the scalar self-energy by introducing a rather 
strong tensor interaction in the model but still retaining 
a reasonable description of other properties. 

Considering the above observations it is natural to ex- 
tend the standard relativistic density functional in a way 
that energy or momentum dependent self-energies ap- 
pear in the Kohn-Sham equations. This will be a phc- 
nomcnological approach introducing new parameters in 
the model that have to be fitted to properties of nuclear 
matter and atomic nuclei. In Ref. couplings of the 
meson fields to derivatives of the nucleon field were in- 
troduced in the Lagrangian density with the desired re- 
sult. In this derivative coupling (DC) model additional 
densities in the functional appear that were not consid- 
ered before. It was possible to describe the experimen- 
tally observed energy dependence of the optical potential 
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by adjusting the relevant coupling constants but still re- 
taining reasonable nuclear matter properties. There are, 
however, some deficiencies of the approach. The density 
dependence of the momentum dependence was different 
for scalar and vector self-energies leading to problems 
in extrapolations of the model to high densities. Since 
the relativistic energy of the nucleon contains the rest 
mass rather large contributions of the energy-dependent 
part of the self-energy are required to show a sizable ef- 
fect. Simultaneously, a large modification of the standard 
energy-independent part of the self-energies was required. 
An attempt to fit the parameters of the model to prop- 
erties of finite nuclei proved to be very difficult. 

Realizing the problems, this paper introduces a mod- 
ification of the model in Ref. 53] combining the virtues 
of the density-dependent approach for the coupling con- 
stants with the momentum-dependent self-energies. In 
Sec. m the basic Lagrangian density of the model with 
density-dependent and derivative couplings (D'^C) is in- 
troduced and the field equations for nucleons, mesons and 
the photon are derived. The relevant equations for nu- 
cler matter and finite nuclei are presented in Sec. lIIII The 
parametrization of the coupling functions and the fit of 
the parameters to properties of finite nuclei are described 
in Sec. IIVI Results of the model for nuclear matter and 
finite nuclei are discussed in Sec. with a comparision 
to conventional RMF models. Finally, in Sec. IVII con- 
clusions and an outlook complete the paper. 



II. LAGRANGIAN DENSITY AND FIELD 
EQUATIONS 

In the present approach the (symmetrized) Lagrangian 
density assumes the form 
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'^PTf,iD^'i: + {iD^'^P)T^^ - ^rM^^ + (10) 



with the nucleon field -0, the covariant derivative 



and the mass operator 

M* =m-T„a -T/,f -5 . 



'-A,. 



(11) 



(12) 



The meson fields are denoted by tr, uj^, 5, and the pho- 
ton field by A^. For completeness all four combinations 
from the alternatives scalar- vector and isoscalar-isovector 
are included for the mesons. The quantities F^, F^, 
Fp and F-,, specify the coupling strenght of the mesons 
and of the photon, respectively, to the nucleon. The el- 
ements of the vector r are the isospin matrices. The 
contribution 
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1 



is the Lagrangian density for the free mesons with masses 
Wcri Wojj and the photon with the field tensors 



F, 



d^Au — diiA^ 
Hpi, = df_,p^ - dyp^ . 



(14) 
(15) 
(16) 



In standard relativistic models the quantities F^ and F 
in Eq. H1U|) are the Dirac matrices 7^ and the unit matrix, 
respectively. In the present approach, they are given by 



F = 



1 + l^u^Y^"" - UpZ^ 



with the quantities 



rn 



4 ^ 



(17) 
(18) 



(19) 
(20) 



that depend on the isoscalar a and w meson fields. An 
extension of the approach to isovector fields is obvious. 
The particular dependence on the meson fields is required 
by the Lorentz structure in order to generate a linear 
energy dependence of the scalar and vector self-energies 
(see below). Furthermore, for both Y^'^ and Z^^ a form 
quadratic in the meson fields was chosen so that the den- 
sity dependence of the energy dependence is similar for 
both self-energies in contrast to the DC model ,52,]. Fy 
and Fg represent the two additional couplings of the D'^C 
model that can depend on the nucleon fields similar as 
the quantities Fcr, T^, Tg, and Fp in the minimal coupling 
of the nucleon field to the meson fields. In Eqs. H17() and 
(fH^ the metric tensor g^i, = diag(l, —1, —1, —1) and the 
four velocity = j^./ Qv depending on the vector current 
density j^j, — V'TmV' with the vector density Qv = \J jp.]^ 
appear. Combining the usual minimal meson-nucleon 
couplings with the action of the quantities F^ and F one 
finds that terms with third powers in the meson fields 
appear in the Lagrangian like in the non-linear models 
with meson self interactions. But in the present approach 
there are no free parameters for these individual terms 
and they are all accompanied by nuclear fields. Sim- 
ilarly to the standard relativistic models with density- 
dependent couplings we assume that the couplings Fo-, 
Ti^, F5, and Fp as well as Fy and F5 depend on the 
vector density g^- In principle, one can also imagine a 
dependence on other densities, e.g., the scalar density 
Qs = tpi^- The explicit choice of the functional form will 
be discussed in Sec. IIVI 

The Dirac equation, i.e. the Kohn-Sham equation, for 
the nucleons in the D'^C model 



7'^ (z9^ - Sp) 7^ - (m - E) V = 



(21) 



5 



has the standard form of relativistic models with the vec- 
tor self-energy 



R 



and the scalar self-energy 



where 



I + T3 



A.. 



(22) 



(23) 



(24) 



and 



s = r^a + Tsf-S--d''Zn 



(25) 



The self-energies in the D^C model are differential opera- 
tors that act on the nucleon field tp. The momentum de- 
pendence enters through the contributions proportional 
to (iD^ — M*u^). The time component of this expres- 
sion mimics the kinetic energy E — m when applied to a 
plane wave in the absence of meson fields. This is in con- 
trast to the DC model where only the iD^ term appears 
p. 

The vector self-energy ()22f) contains the rearrangement 
contribution 



Kuj^J^ + r'p^ ■ > - r>F, - r^-^". Ps (26) 



Qv 



with derivatives T[ 
densities 



dVi/dg^ {i = a, uj, 6, p, V, S) and 



(27) 
(28) 



that replace the usual vector and scalar densities. Addi- 
tionally, the densities 



hi' 



(29) 

(30) 

(31) 
(32) 



appear in the rearrangement contribution 12()|) . ft is 
slightly more complicated than in the standard density- 
dependent (DD) models due to the occurence of the terms 
with y^jy and Z^. ft simplifies considerably for stationary 
systems. The kinetic densities (I29f) and (|3UII are different 



from the derivatives i-D^j^ and iDf^Qs of the vector cur- 
rent and scalar density Qg. The former do not vanish 
in homogeneous nuclear matter and give a finite contri- 
bution. 

From the Dirac equation (|21l) the continuity equations 



a^J^ = and dfJ^ 







(33) 



are derived. Correspondingly, Jq E^nd Jo are the con- 
served isoscalar and isovector baryon densities instead of 
jo = Qv and jb = in standard RMF models. 

In the D'^C model the field equations of the mesons are 
obtained as 



(34) 
(35) 



df,d^'6 + mi5 = TsPs (36) 



with 



and 



(37) 



(38) 



(39) 



The source terms are given by simple products of the 
density-dependent coupling functions To-, T^, F^, Fp, and 
the densities (|27f) and (|28|) . For Fg there is a coupling 
of the equations for the a and lo field. Explicit terms with 
Fy appear only in the equation for the uj field due to the 
factor C^u- The field equation of the photon 



V T 



(40) 



has the usual form with the conserved charge current 
density J^f_, = [J^ + (J^)3]/2. 



III. 



STATIONARY SYSTEMS 



The equations of motion for the nucleons, mesons and 
the photon field simplify considerably if the nuclear sys- 
tem possesses certain symmetries. E.g., in stationary 
systems the meson fields are independent of time. For 
nuclear matter and finite nuclei, it suffices to consider 
only the time-like component of all four-vectors and the 
third component of the isospin vectors. The conserved 
baryon density 



g = Jq= jo (1 + Yqo) - PsZo 



(41) 



depends on the usual vector density jo — Qv and the stan- 
dard scalar density Qs- Similarly, the generalized scalar 
density is given by the combination 



Ps ^ Qs {1 - Zo) + joYoi 



(42) 
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with the quantities 



Yr 



00 



2 2 

nrr- 



and 



-uJocr . 



(43) 



(44) 



Corresponding equations hold for the isovector densities 
g — Jo and Ps- The self-energies in the Dirac equation 
can be written as 



and 



with 



and 



where 



and 



Vb = Wo + YoQ {vQ+m- s) + 
S = s + Zq {vf) + m — s) 



vo = T^ujQ + Tpf ■ po + T. 



1+T3 



^0 



s = r^a + Tsf ■ S 
with the rearrangement contribution 



(45) 
(46) 
(47) 
(48) 
(49) 
(50) 



yli _ 



r>o J° + T'po ■ - KctPs - r'J ■ Ps (51) 



3^ 



r' 



The partial derivative id^ in H45|l and (|46|l gives the 
single-particle energy E when applied to a single-particle 
state in nuclear matter or finite nuclei leading to an ex- 
plicit energy dependence of the self-energies. Note that 
the potentials wq, s, Vq, 5*, and the self-energies Sq and 
S are generally different for protons and neutrons. 



A. Nuclear matter 

Solutions of the Dirac equation are given by the plane- 
wave states 



'/'(P, cr, t) = u(p, cr, r) exp (-ip^a;^) 



(52) 



for a nucleon with four momentum — {E,p), energy 
-E > 0. The positive-energy four-spinor is denoted by 
u(p, (7, r) with spin and isospin quantum numbers a and 
r, respectively. The energy Er of a proton (r = 1) or 



neutron (r = — 1) with momentum p is found in the same 
way as in Ref. |5^j by replacing Xqq with Yooj Yq with 
Zo, and setting W = in the corresponding expressions.. 
It is given by 

Er = (l+Foo)"' (53) 
1 



Vor- 



with the quantities 



ArB + ^Al+p^) 



m — St + BVor 



and 



B = 



l + Yn, 



(54) 



(55) 



The energy density e and the pressure p are calculated 
from the energy-momentum tensor. We find 



mPs 



(56) 



and 



+T^ujf)jQ + Tppo ■ Jo - VadPs - Ts5 ■ Ps 

^[22, 2 ->2 2 2 2 c2 



r 

,-'-[221 2-2 22 2?2 



The densities and gff. are calculated similarly as in 
Ref. They are given by 



(57) 



and 



E _ 2ArBIl + AlB'^Il + II 
~ 7r2(l + yoo)(l- 52)3/2 

M ^ 2ArBIl + Am + 52/J 

^'^^ " ^2(1 +yg„)(l_ 52)3/2 



with the integrals 



^2 — ~ 



-'3 



'p^E^-AllnP^ + ^^' 



Ar 



1 



(58) 

(59) 

(60) 
(61) 
(62) 



c.f. Ref. [5^ . The Fermi momentum p^ is determined 
from = I^'/tt^ by the proton and neutron densities. 

The energy E^ — J A^ + {p!^)^ also depends on the 
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quantity A-r. The standard scalar and the vector den- 
sities can be expressed as 



BIl + Aril 



7r2(l + roo)(l-i?2)i/2 



ArBIl 



7r2(l + roo)(l-i?2)i/2 
We also need the kinetic densities 



and 



with 



''00 



1 + Fc 



00 



1 + ^00 



(63) 
(64) 

(65) 
(66) 



E ^ Ar{l + B^)Il + AlBIl + BIl 

7r2(l + yoo)(l- 52)3/2 ^"'^ 

in the field equations of the mesons. In equation lfF7|) 
for the pressure p a rearrangement contribution appears 
explicitly. It guarantees that the D'^C model is ther- 
modynamically consistent and the Hugcnholtz-van Hove 
theorem jS^J holds as in the case of the DC and DD mod- 
els. 



B. Finite nuclei 

The densities and fields in the ground state of finite nu- 
clei do not vary with time but they depend on the spatial 
coordinates. In the Hartree approximation the state of 
the nucleus is described by a product of single-particle 
states \(t)i) with single-particle energies Si (including the 
rest mass m) that are solutions of the Dirac equation. 
From the nucleon wave functions the various densities 
are calculated by a summation over the contributions 
from the individual nucleons. The meson fields are found 
by solving the corresponding field equations with space- 
dependent source terms with appropriate boundary con- 
ditions. 

The energy of an atomic nucleus in the D^^C model is 
given by 



E = 



+ ^ / [^r- + 2Cjo) <jPs (68) 

+ (r^ + 2T',jo) 5-P,~{T^ + 20o) c^o Jo 
- {Tp + 2T'pjQ) pq- Jq~ T^AqJ^o 



+ / (fr 



(jo^-^^r)(i+^jo)^o 



r' 



with single-particle occupation numbers Wi (in the 
present case without pairing Wi = or 1, respectively). 



Besides the usual rearrangement contributions due to the 
density dependence of the nucleon-meson couplings T^-, 
Ts, r^, and Tp, additional terms appear with the fields 
Yoo and Zq that contain derivatives of the couplings Ty 
and r^. 

Here, we consider only spherical nuclei. The self- 
consistent calculation of the spherical nuclei was per- 
formed in coordinate space in an angular momentum ba- 
sis with a discretization of the radial coordinate and a 
mesh spacing of 0.1 fm using similar procedures as in 
|4l|. A correction of the Coulomb field was introduced 
as in • Because the translational symmetry in the cal- 
culation is broken a cm. correction has to be added to 
the energy H68|) . It has been microscopically calculated 
in the non-relativistic approximation 



Er 



2mA 



(69) 



with the cm. momentum P — X^iLiP** from the single- 
particle wave functions (551] . Nuclear radii are also cor- 
rected for the cm. motion. In case of the charge dis- 
tribution, corrections due to the cm. motion and the 
formfactors of protons and neutrons are considered in 
the calculation of the charge formfactor [2^ . 



IV. 



PARAMETRIZATION 



In order to calculate actual properties of nucler mat- 
ter and atomic nuclei the parameters entering the La- 
grangian have to be specified. For the mass of the nu- 
cleon, the CO and the p meson the conventional values 
TO = 939 MeV, = 783 MeV and = 763 MeV were 
used. The mass of the a meson m„ is treated as a free 
parameter. The S meson is neglected in the present work. 
Except for the electromagnetic coupling constant F^, all 
parameters in the coupling functions F; and the mass 
of the a meson TOo- were obtained from a fit to proper- 
ties of the eight doubly magic spherical nuclei ^^O, 2^0, 
^°Ca, 48ca, 56Ni, looSn, ^^"^Sn, 208pb that include nuclei 
close to and far from the valley of stability. Other nuclei 
were not considered because effects due to nucleon pair- 
ing have to included explicitly in this case adding more 
parameters in the fit. 

There is no unique and generally accepted strategy for 
the fitting procedure 0| . In the present fit of the pa- 
rameters experimental data for binding energies, charge 
radii, diffraction radii, surface thicknesses, and spin-orbit 
splittings were taken into account ^22,, .23] . Furthermore, 
it was required to reproduce the experimental value of 
0.20(4) fm for the neutron skin thickness in ^^^Ph, i.e. 
the difference r„ — rp of the neutron and proton rms radii 
[stI - In total there are 32 experimental data points used 
in the fit. 

For a comparison of the D'^C model with standard 
density-dependent models, a new parametrization for DD 
models was derived under the same conditions in the fit 



8 



as the D^C model with the same set of miclei and experi- 
mental data. In this fit the nuclear incompressibility was 
fixed to K — 240 MeV as in the original model in Ref. 
|4l| where only nuclear binding energies were considered. 
A fit without this constraint leads to unacceptably large 
values of K w 300 MeV. 

The D^^C model with momentum-dependent self- 
energies contains two additional coupling functions, Ty 
and Fs, as compared to usual DD models. They allow to 
fit additional properties of the nuclear system under con- 
sideration. Here, it was required that the optical poten- 
tial in symmetric nuclear matter at saturation density 
assumes the value 50 MeV at a nucleon energy of 1 GeV, 
a value that is typical for Dirac phenomenology 0, ^3 ■ 

The density dependence of the coupling functions is 
written as 



(70) 



with the coupling constants Ti(g,:cf) at a reference density 
Prof and suitable functions fi that depend on the ratio 



Prof 



(71) 



with the vector density p„. Note that the vector density 
in the D'^C model is not identical with the baryon density 
g. In the DD model the reference density is just the 
saturation density Psat of symmetric nuclear matter. The 
reference density in the D^C model is different from psat 
but it corresponds to the vector density g^ determined 
at saturation of symmetric nuclear matter. 

The functional form of fi (x) for the density dependence 
of the a, Lo and p meson is assumed to be the same as 
introduced in Ref. 41] and later used in the parametriza- 
tions DD-MEl (42], DD-ME2 "is^, and PKDD For 
the <7 and w meson it is given by the rational function 



l + c,{x + d,y 



(72) 



with constants a^, 6^, q, and di. In order to reduce the 
number of free parameters it is required that the func- 
tions fa and fi^ obey the conditions /(t(1) = /uj(l) = 1, 
faiO) = /:(0) = 0, and /;'(1) = /^'(l). For the p meson 
a simple exponential law 



(73) 



fp{x) =exp [-Up (a;-l)] 

is assumed with one free parameter Cp as in Ref. 

In the D'^C model the quantities Foo and Zq are re- 
sponsible for the energy dependence of the self-energies. 
Assuming constant couplings Ty and F5 both Yqq and 
Zq show an approximate g'^ dependence at low densities 
g. This suggests to introduce the power law 



as the functional form in order to be able to arbitrarily 
adjust the density dependence of the energy dependence. 



TABLE I: Mass of the a meson rria, parameters of the cou- 
pling functions and reference density g^^i in the DD and D^C 
model. 





UD 




rria [MeV J 


547.204590 


555. 986205 




10.D55257 


ii.(j27474 


aa 


1.371545 


i.o4Do41 


7, 

Oct 


r\ CI A A CvG 

U.d44Uoo 




Co- 


1 no /I c CO 


'7 r\'v\ 'van 
(.uiL ryy 




0.5d7o27 


0.217107 




1 Q 1 ooon 


13. /O0o49 




i. 00000 ( 


2. 004516 


7 


n col ^70 A 
U.ozi 1 z4 


4.563(03 




0.869983 


9.864792 




0.618991 


0.183821 




3.639023 


3.917462 


ap 


0.4987 


0.4220 




0.0 


302.188656 


rs(£'rcf) 


0.0 


-21.632122 


Qrci [fm^''] 


0.148746 


0.128941 



In this work the parameters ay and as are set to unity 
because in this case the second integral in the energy 
(|68|l of an atomic nucleus does not contribute and the 
calculation is simplified. 

In total there are ten free parameters in the D'^C model 
and eight in the DD model. Since several of the pa- 
rameters are strongly correlated f5(| it is convenient to 
perform the fit not directly in all of these parameters 
directly but to use saturation properties of symmetric 
nuclear matter (see subsection lV A|l as independent vari- 
ables in the fit and convert these quantities to the cou- 
pling constants analytically. Considering the conditions 
for the neutron skin thickness in 208p|-, ^^^^^ ^j^g optical 
potential, the minimum of has to be found in an eight- 
(seven-) dimensional parameter space. 



V. RESULTS 

The actual values of the model parameters as deter- 
mined in the fit are given in Table The main differ- 
ence between the DD and the D^C model are the values 
of the parameters a^, 6^, q, and di for i — a^uj in the 
functional form of the density dependence l|72|l . This 
corresponds to a much stronger increase of the a and 
Lu coupling for g ^ in the D^C model as compared 
to the DD paramctrization. The density dependence of 
the p meson, however, is smaller for D'^C than for DD. 
Obviously, the constants Tyigrei) and rs(pref) are non- 
(74) vanishing only in the D'^C model. Note that Fy is more 
than ten times larger than Fg. Correspondingly, the en- 
ergy dependence of the vector self-energy is much larger 
than the energy dependence of the scalar self-energy in 
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TABLE II: Contributions of nuclear properties with assumed 
uncertainty in the fit of the model parameters and total 



80 



property 


uncertainty 


DD 


D^C 


binding energies 


0.2 MeV 


225.2 


168.4 


charge radii 


0.01 fm 


44.9 


32.4 


diffraction radii 


0.01 fm 


74.7 


53.3 


surface thicknesses 


0.005 fm 


15.0 


23.0 


spin-orbit splittings 


0.2 MeV 


95.4 


60.3 


total 




455.3 


337.4 



the present parametrization of the D'^C model. The neg- 
ative sign of r5(£ii.of ) indicates that the scalar self-energy 
rises weakly with the energy of the nucleon. In con- 
trast the vector self-energy decreases much more strongly. 
From Dirac phenomenology 23| one would expect 
that both the scalar and the vector self-energy decrease 
at larger nucleon energies. However, by calculating the 
momentum-dependent contributions to the self-energies 
from the Fock terms considering one-pion exchange (that 
is not included in conventional RMF models), an oppo- 
site energy dependence of the self-energies is expected 
[ssf with an increase of the scalar self-energy and a de- 
crease of the vector self-energy. The result of the D'^C 
model lies between these two cases. 

In Table ^1 the contributions to the total by the 
various nuclear properties in the fit of the parameters 
are given for the DD and D'^C model, respectively. In 
both cases the error in the binding energies contributes 
approximately one half to the total x^- The D^C model 
gives a better description of all properties except for the 
surface thickness with a small increase of the correspond- 
ing partial . The total of the D^C model is about 
27 % smaller than in the DD model. This represents a 
considerable improvement in the description of the ex- 
perimental data. More details for finite nuclei will be 
discussed in subsection IV Bl 



A. Nuclear matter 

The equation of state (EOS) of symmetric nuclear mat- 
ter, i.e. the binding energy per nucleon as a function of 
the density p, is depicted in Figure Q The EOS in the 
D^C model and the DD model are compared to the result 
for the nonlinear NL3 parametrization j32j that is widely 
used in RMF calculations with considerable success. All 
three curves are very similar below a density of approx- 
imately 0.2 fm~^. At high densities there is a noticable 
difference. The parametrization NL3 leads to the stiffest 
EOS and the DD model has the softest EOS with the 
D'^C curve lying in between. 

A more quantitative comparison of the models is pro- 
vided by examining the characteristic constants in the 
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FIG. 1: Equation of state for symmetric nuclear matter in 
various parametrizations. 



expansion of the binding energy per nucleon 
E 



-j = ay 



K K' 

162 



18 



(75) 



L 



K 



18 



in nuclear matter near saturation j59l l63| . The deviation 
of the density from saturation is quantified by 



Q - Qsat 
e = 

Qsat 

and the neutron-proton asymmetry is given by 



Qn 



Qp 



(76) 



(77) 



In symmetric nuclear matter, i.e. (5 = 0, the binding en- 
ergy at saturation oy, the incompressibility K^, and 
the derivative K' of the incompressibility determine the 
form of the EOS. For asymmetric nuclear matter, also 
the symmetry energy J, the derivative L of the symme- 
try energy, and the symmetry incompressibility fTsym are 
important. The values of these quantities for the three 
models NL3, DD, and D^C are given in Table |lll| A de- 
tailed comparison shows that there are major differences 
between the parametrizations. The D'^C model has the 
softest EOS near saturation (i^oo) but at higher densi- 
ties is becomes stiffer than the DD model because of the 
large negative K' . The saturation density is larger than 
in the DD and NL3 models where the latter model also 
has a rather strong binding. 
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TABLE III: Properties of symmetric nuclear matter at satu- 
ration in various models. 





NL3 


DD 


D^C 


0<i^f ffrn "^l 


0.1482 


0.1487 


0.1510 


av [MeV] 


-16.240 


-16.021 


-15.981 


Kr^ fMeVl 

*■ OO L-LViC V J 


271.5 


240.0 


232.5 


K' [MeV] 


-203.0 


-134.6 


-716.8 


J [MeV] 


37.4 


31.6 


31.9 


L [MeV] 


118.5 


56.0 


59.3 


Ti-sym [MeV] 


100.9 


-95.3 


-74.7 


m* /m 


0.596 


0.565 


0.541 


meff/m 


0.655 


0.628 


0.710 


K:cn(l GeV) [MeV] 


282.2 


310.0 


50.0 




0.1 0.2 0.3 0.4 0.5 

density p [fm ] 



FIG. 2: Symmetry energy J in symmetric nuclear matter as a 
function of the nucleon density g in various parametrizations. 



The symmetry energy J and the derivative L are con- 
siderably smaller in the D'^C and the DD model than 
in the NL3 parametrization. This is due to the density 
dependence of the p meson coupling. The actual values 
are determined by the fit of the neutron skin thickness 
in ^"^Pb. The differences in L corresponds to a rather 
strong deviation in the density dependence of the sym- 
metry energy that is shown in Figure [5] In the NL3 
model the symmetry energy rises almost linearly with 
the density. In contrast, the DD and D'^C model exhibit 
a considerable flattening. This behavior is also reflected 
in the EOS for neutron matter as presented in Figure 
El Here, the NL3 model shows a very stiff EOS. The 
shape of the neutron EOS in this model also differs at 
low densities from the results in the DD and D^C model 




0.1 0.2 0.3 0.4 0.5 

density p [fm ] 



FIG. 3: Equation of state for neutron matter in various 
parametrizations. 



that display very similar results for neutron matter. It is 
worth noting that the symmetry energy and the neutron 
EOS in the three models are rather similar at a density 
near 0.1 fm~'^. This value is approximately the neutron 
density in heavy atomic nuclei. It would be rewarding to 
study the effects of the different equations of state on the 
properties of neutron stars but this is beyond the scope 
of the present work. 

Besides the parameters characterizing the EOS, also 
the Dirac mass m* and the Landau mass rricg (at the 
Fermi momentum) are of great interest. The values of 
these quantities for saturated symmetric nuclear matter 
are also given in TableHni The D^C model has the small- 
est value for m* even below that of the DD model but 
the Landau mass of 0.710 nucleon masses is the largest 
of the three models. The of the fit is actually not very 
sensitive to a change in the Landau mass. In principle 
it is possible to adjust rrieS to even larger values at the 
Fermi surface in saturated symmetric nuclear matter in 
the D^C model. 

The density dependence of m* and nics in symmet- 
ric nuclear matter is shown in Figure 01 The relativis- 
tic effective mass drops strongly with increasing density. 
At low densities the D^^C model is similar to the DD 
model whereas at high densities it follows more closely 
the NL3 model. At low densities the Landau mass de- 
creases like the relativistic effective mass with increas- 
ing density. However, at higher densities an increase is 
observed that is easily explained considering the depen- 
dence on the Fermi momentum in equations © and Q . 
In the D'^C model this increase is very pronounced and 
the absolute values of Wcff are much larger than in the 
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FIG. 4: Relativistic effective mass m* (a) and effective Lan- 
dau mass moft at the Fermi momentum (b) in units of the 
nucleon mass m in symmetric nuclear matter as a function of 
the density in various parametrizations. 



FIG. 5; Schrodinger-equivalent central optical potential Veen 
in symmetric nuclear matter as a function of the nucleon ki- 
netic energy iSkin — E — m in various parametrizations and 
from Dirac phenomenology (DP) of elastic proton-nucleus 
scattering [iali^. 



other parametrizations that are representative for stan- 
dard RMF models. Correspondingly, the density of states 
at the Fermi surface in the D'^C model is considerably 
higher than that in other RMF models. From the expe- 
rience in the fitting of the parameters it is found that the 
requirement of a larger Landau mass at saturation than 
in the present D'^C model leads to an even softer equation 
of state for symmetric nuclear matter if one also demands 
an optical potential of 50 MeV at 1 GeV kinetic nucleon 
energy. Thus, this last condition should be relaxed in a 
more extensive parameter fit. 

Another major advantage of the D^^C model is the 
possibility for a very reasonable description of the 
Schrodinger-equivalent central opical potential as de- 
picted in Figure m At low kinetic energies of the nucleon 
it rises nearly linearly as the empirical optical potential of 
nuclear matter extracted from Dirac phenomenology for 
elastic proton- nucleus scattering 46, 47] . At higher nu- 
cleon energies it shows a similar saturation with reason- 
able absolute values. In standard RMF models like NL3 
and DD without momentum-dependent self-energies the 
optical potential rises linearly with the nucleon energy. It 
approaches unrealistically high values around 300 MeV 
at 1 GeV kinetic energy (see also Table 1111}) . Because the 
optical potential in the D^C model is a quadratic func- 
tion of the nucleon energy, it will ultimately decrease at 
higher energies und become unrealistic. However, this 
behavior is not relevant for most applications, e.g. finite 
nuclei or high-density nuclear matter, because the rele- 
vant Fermi momenta and corresponding nucleon energies 



TABLE IV: Total binding energies (in MeV) of the eight fit 
nuclei in the experiment [6j| and in the models DD and D"^C. 



Nucleus 


Exp. 


DD 


D=^C 




-127.619 


-128.064 


-127.040 


24 Q 


-168.382 


-169.114 


-169.178 


*°Ca 


-342.052 


-342.505 


-342.581 


*«Ca 


-415.990 


-414.876 


-415.047 


56 Ni 


-483.992 


-481.924 


-482.486 


"'«Sn 


-824.794 


-826.255 


-826.079 


"^Sn 


-1102.851 


-1103.359 


-1103.478 


208pb 


-1636.430 


-1636.030 


-1635.893 



are still small enough in these cases. 



B. Finite nuclei 

From Table ITU it was already seen that the D'^C model 
improves the description of various properties of finite nu- 
clei as compared to the DD model demanding the same 
conditions in the fit. In Table ITVl the experimental bind- 
ing energies of the nuclei that were considered in the 
fit are compared with the results of the DD and the D'^C 
model. In general, there is a good reproduction of the 
experimental data. The absolute root-mean square devi- 
ation from the experiment is 1.06 MeV and 0.92 MeV, 
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FIG. 6: Charge radius (a), diffraction radius (b), and surface 
thickness (c) of the fit nuclei in the DD model (blue diamonds) 
and in the D'^C model (red squares) compared with data from 
the experiment (solid black dots) |65i.i66^| . 



respectively. These values are still larger than corre- 
sponding numbers from dedicated fits to masses with rms 
deviations in the order of a few hundred keV, see, e.g., 
[60, 61, _6^]. For a fair comparison one has to bear in 
mind that the RMF parametrizations in this work are 
not just mass models. They are constructed in order to 
describe a large number of different properties of nuclei 
and nuclear matter. Furthermore, a valid comparison is 
only possible when a larger set of (also deformed) nuclei 
is taken into account. However, this requires a correction 
for the rotational energy and a reasonable description of 
pairing effects that is left for future investigations. 

The characteristic parameters of the nuclear shape and 
of the charge form factor, respectively, i.e. the charge 
radius, the diffraction radius, and the surface thickness 
[2^ |23 |. are compared to experimental data [6^ 
in Figure |B| Both the DD model and the D'^C model 
agree very well with the available experimental data. At 
the same time the difference between the models is very 
small. They show the same trend for the radii and for 
the surface thickness. The largest deviation from exper- 
iment is found for the surface thickness of ^^O, a light 
nucleus that is not expected to be described very well in 
a mean-field model. 

The last quantity that was considered in the fit is the 
spin-orbit splitting of proton and neutron single-particle 



TABLE V: Spin-orbit splitting (in MeV) of neutron {u) and 
proton (tt) levels in the experiment and in the models DD and 
D^C. 



Nucleus 


State 


Exp. 


DD 


D^C 


16 Q 


i^Op 


6.18 


6.761 


6.278 




TrOp 


6.32 


6.695 


6.218 


4«Ca 


uOi 


8.39 


7.961 


7.623 




vlp 


2.03 


1.512 


1.585 




vQi 


7.16 


8.568 


8.369 




ulp 


1.12 


1.445 


1.335 


Sn 


TTld 


1.74 


1.996 


1.786 




TfOg 


6.09 


6.563 


6.171 


2'«Pb 


i/2p 


0.90 


0.917 


0.879 




TTld 


1.33 


1.834 


1.615 




TrOh 


5.56 


6.063 


5.642 



energies of levels with the same principal quantum num- 
ber and the same orbital angular momentum. In Table Ivl 
the results in the DD model and the D'^C model are com- 
pared to experimental data for levels close to the Fermi 
surface. It is found that the spin-orbit splittings of the 
D^C model improves the description in most cases when 
compared to the DD model. This is possible because 
the relativistic effective mass m* in the former model is 
even smaller than in the latter model. At the same time, 
the Landau mass TOoff in the D'^C model is considerably 
larger than in the DD model. Correspondingly, the level 
density of the D'^C model is higher. This compression 
of the spectrum is easily seen when the distribution of 
single-particle levels is compared for the three heaviest 
nuclei in the fit in Figured Levels below the Fermi en- 
ergy are shifted to higher energies and levels above the 
Fermi energy become more strongly bound. 



VI. 



CONCLUSIONS AND OUTLOOK 



In this work the conventional relativistic mean-field 
model with density-dependent nucleon-meson couplings 
was extended by introducing a particular form of cou- 
plings between the isoscalar meson fields and derivatives 
of the nucleon fields. This approach leads to a linear mo- 
mentum dependence of the scalar and vector self-energies 
in the Dirac equation for the nucleon. The parameters 
of the model were determined by a fit to properties of fi- 
nite nuclei. It was possible to improve the description of 
binding energies, nuclear shapes and spin-orbit splittings 
of single-particle levels. The characteristic parameters of 
nuclear matter in the D'^C model are shifted closer to 
the values of non-relativistic Skyrme Hartree-Fock mod- 
els [63. Most noticable was the increase of the effective 
Landau mass and, correspondingly, the level density at 
the Fermi surface as compared to standard RMF mod- 
els. At the same time, the momentum dependence of the 
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FIG. 7: Neutron (i/) and proton (n) single particle energies 
for (a) ^°"Sn, (b) ^^^Sn, and (c) ^^^Pb in the DD (blue lines, 
left) and D^C (red lines, right) model. The position of the 
shell closure is denoted by open circles. 



self-energies leads to a Schrodinger-equivalent optical po- 
tential in symmetric nuclear matter that follows closely 
the empirical data whereas standard RMF models fail. 

The introduction of couplings of the meson fields to 
derivaties of the nucleon field is a purely phcnomenolog- 
ical approach. The goal was to investige possible exten- 
sions of the RMF model and their effect on various quan- 
tities in nuclear matter and finite nuclei. The origin or 
the underlying mechanism of the additional couplings is 
not really relevant in the present context. They are seen 
only as an effective way to introduce some desired effects 
in the spirit of density-functional theory. This approach 
can be compared with the medium dependence that can 
be parametrized in relativistic mean field models by self- 
interactions of the meson fields in non-linear models, by 
density dependent couplings, or even by density depen- 
dent meson masses. These approaches are more or less 



equivalent when the final results for properties of nuclei 
are compared, however the underlying mechanisms are 
quite different. 

In principle, of course, one could go beyond the mean 
field model and discuss in a systematic diagrammatic ex- 
pansion modifications of the self energies. But this was 
not the topic of the work. In the present paper, I wanted 
to stay on a purely phencmcnological level in a simple 
self-consistent approach that has all the necessary fea- 
tures that make it possible to apply the model in practical 
calculations with reasonable effort to nuclear matter and 
finite nuclei. The results show that the approach works 
surprisingly well. The modification of the Lagrangian 
density can be discussed and treated without reference 
to more fundamental mechanisms that could generate the 
corresponding terms. 

The density dependence of the momentum dependence 
is determined in the present model by the choice of the 
parameters ay and as in Eq. H74(l . Other values than 
av — as — I SiS in the present model should be consid- 
ered. A fit of the model parameters by fixing the Landau 
mass at the Fermi momentum but varying the absolute 
value of the Schrodinger-equivalent optical potential is 
conceivable. The parameter space can be explored more 
thoroughly allowing for a simultaneous variation of the 
optical potential and of the Landau mass. The study of 
deformed nuclei is a further test of the model. However, 
it will require a larger computational effort and the con- 
sideration of pairing effects. It will give more insight into 
the quality of the description. Applications to neutron 
stars and to simulations of heavy-ion collisions are also 
possible. Extensions of the model to introduce an addi- 
tional isospin dependence of the momentum dependence 
are obvious. 
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